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A DEVICE INDEPENDENT COMPUTER PLOTTING SYSTEM 



UDAY G. GUJAR 

Computing Centre and 
School of Computer Science 
University of New Brunswick 
Fredericton, N.B., Canada 
E3B 5A3 



1 . 0 Introduction : 

It is an established fact that the 
human mind comprehends and understands 
pictorial information very easily and 
rapidly. With this experience and the 
advent of computers r it is little sur- 
prise that the computers have been and 
are being used more and more to produce 
graphical output. 

Almost every computing centre, may 
it be small, medium or large, has some 
sort of plotting hardware and the soft- 
ware to support it. More often than 
not, the software is device dependent. 
As a result, while switching from one 
plotting device to another, the user is 
burdened with several programming 
changes and forced to live with the 
associated delays and inconveniences. 
Further, addition of a new plotting 
device is a traumatic experience both 
for the users as well as the system 
programmers who implement the soft- 
ware for such a device. 

This paper describes a computer 
plotting system which is completely 
device independent. The user can 
switch from one plotting device to an- 
other without any programming changes. 
The interface for adding new plotting 
devices is formalized and discussed. The 
addition of a new plotting device is com- 
pletely transparent to the user; in fact, 
the old programs may be used to produce 
plots on the new device without any pro- 
gramming changes whatsoever. 



The system was designed and implemen- 
ted by the author in late 1971 and has 
been in use- ever since (see Ref . 1) . 
Several new plotting devices, real as well 
as pseudo, have been added successfully 
without affecting the users of the system. 



2.0 Devices and Modes : 

2 . 1 Available Plotting Devices : 

Plots can currently be obtained on 
the following plotting devices: 

a. 1403 line printers (set to print m 
characters per inch and n lines per 
inch where m and n may have any 
values) . 

b. 611 Tektronix storage oscilloscope. 

c. 1627 Calcomp plotters. 

d. 5100 Gould electrostatic plotter. 

It is interesting to note the diverse 
functional characteristics of the above 
plotting devices. The 611 storage oscil- 
loscope is an analog device while the 
1627 Calcomp plotter is a digital incre- 
mental plotter. However, both of these 
devices allow positive and negative move- 
ments in both the X and Y directions. The 
5100 electrostatic plotter and line 
printers have different characteristics 
compared to the 611 or the 1627 in that 
they allow paper movement only in one 
direction. As a result, the programming 
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technique involved is quite different for 
these later devices. 



Step A : Executed by the user. 



The interface between the plotting 
devices and the main computer (currently 
IBM 370/158) are also of diversified 
nature. For example, one 162 7 plotter is 
connected through an IBM 1827 data control 
unit (which is basically a digital to ana- 
log and analog to digital converter) and 
is situated in the same room as the main 
computer. However, another 1627 plotter 
is connected via a HASP remote work sta- 
tion and is situated 130 miles away from 
the main computer. 

Various other plotting devices, for 
example Hewlett Packard X-Y recorders, 
have also been connected to the system - 
though these are not made available to the 
users because of the non-availability of 
these devices at the author's installation 
and such reasons. 

2.2 Plotting Modes : 

Since the plotting usually requires 
some manual initial operations involving 
hardware setup like switching the device 
on, ensuring plotting pen conditions, 
readying the interface, etc., two modes of 
plotting are provided, namely off-line and 
on-line. 



2.2.1 Off -Line Mode of Plotting : 

In this mode, the information 
required to produce a plot is written on a 
storage medium (e.g. a disk or a tape file) 
by the plotting package. At some later 
time, a control program is used to retrieve 
this information to create the plots on any 
of the available plotting devices. This 
process is depicted in Fig. 1. 

Off-line mode of plotting can be 
thought of as a pseudo plotting device. 
Several types of such pseudo plotting de- 
vices are available. These are: 

a. MASSPLOT - This is based on a direct 

access data set which can hold 
several plots of several users 
at the same time. The data set 
management is handled by the 
plotting system. 

b. USERPLOT - This enables the user to 

use his own sequential data set, 
either on a disk or tape, to 
store a number of his plots. 

c. HASPPLOT - This uses the HASP spool 

pack to store the plots and is de- 
pendent upon a local modification 
made to HASP (see Ref. 2). 

There is a control program which cor- 
responds to each of these pseudo devices. 
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Step B : Executed by the Computing Centre 
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Fig. 1: Off -Line Mode of Plotting 



2,2.2 On-Line Mode of Plotting : 

In on-line mode, the user pro- 
gram drives the plotting device directly 
in real time as can be seen from Fig. 2. 
Since about the only justifiable reason 
for doing on-line plotting is to permit 
the programmer to interact directly with 
the plotting programme, the user must 
arrange for his job to be run at a special 
time. He or she would normally supervise 
the plotting personally. 
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Pig. 2: On-Line Mode of Plotting 



2.3 "Simulated" Device ; 

it is recognized that several instal- 
lations support the calcomp basic software. 
To facilitate the transportability of pro- 
grams from such installations, an inter- 
face has been written which simulates 
basic calcomp software. Using this inter- 
face, which is invoked by a single JCL 
card, the programs that use basic calcomp 
software can obtain plots on any of the 
available plotting devices, either dir- 
ectly or through the pseudo plotting de- 
vices, without any changes whatsoever. 

An interesting corollary of the above 
is to simulate this plotting system under 
discussion to generate calcomp basic soft- 
ware calls. Once this is done, the pro- 
grams written at the University of New 
Brunswick, using our plotting system, 
could be sent to any other installation 
supporting calcomp basic software. It 
is the author's belief that this can be 
done rather easily. 



3.0 User's Overview : 

The system is based on several sub- 
routines accessible from FORTRAN IV. 
Effort is made in this section to give an 
insight into this plotting package. A 
brief description of all the device inde- 
pendent subroutines is given. 



3.1 An Example : 

Consider the problem of plotting the 
equation: 

y = sin(x) + sin(2x) for 0sx*2it 



The following program, which is self ex- 
■ planatory, produced the plots given in 
Figs. 3 and 4. 

PI2 = 2.*3. 14159 

C PLOT A HEADING 

C {DEFAULT: 1" = 1 USER UNI T) 

CALL THICMO.iS.ZSt'EXAMPLEn', 

* 0. f .5.6t4) 

C FIX PHYSICAL SIZE TO V BY 5" 

CALL AREAU.iS.) 
C DEFINE USER UNITS 

CALL SETPLT(0.,-2.,PI2,2.) 
C DRAW A GRID 

CALL GRID(0.,-2.,PI2,2.,-*.i-4.l 
C PLACE THE PEN AT INITIAL POINT 

CALL NOWPLTtOfO. ,0.) 
C USE * FOR PLOTTING THE CURVE 

CALL PRNTCHt ) 
C PLOT THE GIVEN EQUATION 

DX=PI2/60. 

DO 1 1=1,60 
X=DX*I 

Y=SIN(X)+SIN(2.*X) 
1 CALL NOWPLT(lrX.Y) 

C IDENTIFY THE CURVE 

CALL CHRPRT(1.4tl.5t 

* l iLY»SINm+SIN(2X)n l »0. i .15) 
C TERMINATE THE PLOT 

CALL ENDPLT 

STOP 

END 

Note that n in the above statements repre- 
sents a multiple punch of 0, 2 and 8. 

3.2 Device Selection : 

As far as the user is concerned, he can 
switch from one plotting device to another 
by means of a simple change in his job con- 
trol language (JCL) . This often involves 
changing one JCL card. 

All the common plot routines and the 
routines for generating printer plots are 
stored in the data set UNB1 • FORTLIB . The 
routines for various pseudo plotting de- 
vices are stored in various different 
libraries, e.g. the routines. required for 
MASSPLOT "device" are stored in 
UNBl-PLOT'MASSPLOT and those that are re- 
quired for USERPLOT "device" are stored in 
UNBl-PLOT'USERPLOT. The routines required 
for on-line plotting are stored in the 
data set UNBl • PLOT • REAL PLOT . 

The proper device is selected by 
painting the SYSLIB DD statement of the 
linkage editor or the loader step to the 
corresponding data set and UNB1 -FORTLIB in 
that order. 

Since most of the users are not very 
experienced with JCL, several procedures 
are written for convenience. Thus, for 
printer plots one may use the following 
JCL: 
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//. . . JOB. . . 

// EXEC FORTGCLG 

//FORT-SYS1N DD * 

... FORTRAN PROGRAM 



For using MASSPLOT, one simply changes the 

// EXEC FORTGCLG 
card to 

// EXEC MASSPLOT 

Similarly for using USERPLOT, above card 
is changed to 

// EXEC USERPLOT 



3.3 Optional Arguments : 

Several routines in this package have 
a variable number or arguments. Many of 
these are completely written in FORTRAN 
(see Appendix B) using the technique de- 
veloped in Ref. 3. Very briefly, from the 
user's point of view, the user may specify 
as many of these optional arguments as de- 
sired. The arguments which are omitted get 
the predefined default values. Arguments 
can only be omitted from the right hand end 
of the list. For example, consider the 
following routine: 

CALL ELL I PS ( XC , YC , RMAJOR [, RM] NOR , THETA , FROHTH , TOTH] ) 

where last four arguments are optional. The 
following calls are valid: 



and a data definition card (DD card) is 
provided to give the attributes of the 
user's data set which might be: 

//G0-FT23F001 DD DI SP=SHR , DSN=MYFILE 



CALL ELLIPS(XC,YCtRMAJOR> 

CALL ELLIPS(XCiYC,RMAJOR,RMINOR) 

CALL ELL I PS { XC , YC » RMAJOR , RM I NOR , THE T A ) 

CALL ELLT PS { XC , YC , RMAJOR ,RM INOR , THETA r FROMTH ) 

CALL ELLI PS t XC t YC t RMAJOR , RM INOR i THETA f FROMTH , TOTH ) 



• • • • • 
• • » • 

• • • • • 



.* 

* 

* 
* 



* • 

• • « • 
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* 
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♦ 

* 4 
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« 
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EXAMPLE 




Fig. 3: Plot Produced on a Line Printer 



Fig. 4i Plot Produced on the 5100 Electrostatic 
Plotter 
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Note that if a special value for THETA is 
required, a value for KM I NOR must be spe- 
cified even if the default value of RMINOR 
is suitable. Similarly, to be able to 
specify TOTH, the values for RMINOR, THETA 
and FROMTH must be supplied in that order. 



3.4 Common Plot Routines and Their 
Functions ; 

All the common plot routines have 
been categorized as far as possible and a 
brief description of their usage is given 
in this section. See Fig. 5 for a pictor- 
ial overview. The user may use any or all 
of these routines for plotting on any de- 
vice. If any particular routine is not 
applicable, by virtue of its peculiarity, 
to some plotting device a dummy routine is 
automatically supplied. This eliminates 
the necessity of making changes, which 
otherwise might have been required , to the 
deck while switching from one plotting de- 
vice to another. See Appendix A for the 
summary of calling sequences of all the 
routines and Ref . 1 for the complete de- 
tails. Appendix B contains such details 
as the memory required for each routine, 
the language each routine is written in, 
entry points and the other routines called. 



3.4.1 Initialization Routines : 

These routines should be called, 
if at all, before any other plot routine. 
However, one may call the subroutines AREA 
and SETPLT several times to achieve some 
fascinating results. The following is a 
list of the initialization routines: 

a. DEVICE: Enables the user to choose be- 

tween the plotting devices (e.g. 611, 
1627 or 5100) during the on-line 
mode of plotting. 

b. AREA: Defines the physical size of the 

plot. 

c. SETPLT: Enables the user to define his 

own units for plotting rather than 
forcing some arbitrary units. 

d. PLOT ID: Enables the user to supply an 

80 character identification label to 
his plot when using the off-line 
mode of plotting. This identifica- 
tion label, which may contain some 
special instructions, is conveyed to 
the plotter operator by the plotting 
system just before the initiation of 
the start of the actual plot. 



Initialization Routines 

DEVICE 

AREA 

SETPLT 

PLOTID 

CHLNIK 

PRNTCH 



Termination Routine 
ENDPLT 



Used only 
for off- 
line 

plotting 



Control Programs 




MASSPLOT 




USERPLOT 




etc . 





Basic Plot Routines 



PLTRTN 
LCTRTK 
SFTRTN 





Character Plotting 

CHARS 
SPSYNB 
THICK 
CHRPRT 



1 



"Nonplotting" Plotting 
Utilities 

LOCATE 
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Plotting Utilities 



Curve Plotting 
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PLTPLT 



Geometric Figures 
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Miscellaneous 
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Fig. 5: Overview of the Plotting Package 
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e. CHLNIN: Defines the number of charac- 

ters per inch and the number of 
lines per inch for the line printer. 

f. PRNTCH: Enables the user to select a 

print character for plotting on a 
line printer; this character may be 
changed as often as required as the 
plot progresses. 



3.4.2 Termination Routine : 

There is one termination routine, 
namely ENDPLT. All the plotting jobs must 
call this subroutine ENDPLT at least once 
to terminate the current plot frame and 
thus force completion of the plot on the 
desired plotting device. 



3.4.3 The User's Pencil : 

The following routines are pro- 
vided to move the plotting pencil: 

a. NOWPLT: Moves the plotting pencil 

from its current position to a 
specified location in a straight 
line either in up or down position. 
This permits the user to draw a 
point or a straight line or to 
position the plotting pencil. 

b. PLTPLT: Enables the user to make 

several sets of plotting pencil 
movements, as specified above in 
"NOWPLT" , in one call. 

c. SPSYMB: See Section 3.4.4. 

d. DSHLNS: See Section 3.4.6. 



3.4.4 Character Plotting : 

A variable width character gen- 
erator has been written which enables the 
user to plot a string of characters con- 
taining all the EBCDIC graphic symbols 
which include upper and lower case letters, 
several special symbols, digits etc. In 
addition, several useful plotting symbols 
are also available. Special control 
characters are provided for writing both 
upper and lower case letters, superscripts, 
subscripts and for starting new lines. 
The following is a list of these routines: 

a. CHARS: Plots a string of characters 

with any specified height at any 
specified angle. 

b. SPSYMB: Plots one of fifteen special 

symbols with specified height at 
any angle. The user may optionally 
join all or any of these points to 
obtain one or more curves . 

c. THICK: Enables the user to plot char- 

acters with specified thickness, 
height and inclination. 



d. CHRPRT: PRINTS, rather than plots, a 
string of characters on a line 
printer. For other devices, this 
routine simply calls 'CHARS' to PLOT 
the characters. 



3.4.5 "Nonplotting" Plotting 
Utilities : ■ 

There are two routines which do 
not produce plots but are useful: 

a. LOCATE: Returns the current position 

of the plotting pencil. 

b. SHIFT: Enables the user to shift, by a 

specified amount, all future plot- 
ting with respect to plotting done 
prior to the call to SHIFT. 



3.4.6 Plotting Utilities : 

It is anticipated that some of 
the sections of plots required by users 
will be recurring in nature. Some examples 
are to draw axes, grids, circles, ellipses, 
to obtain curves joining several points, 
etc. Several routines have been written to 
fulfill this anticipated need, to avoid 
possible duplication of efforts and to pro- 
vide the user with an easy and convenient 
tool. 

A. Curve Plotting : 

a. LNSPLT: Draws straight lines be- 

tween several specified points in 
one call. 

b. DSHLNS: Joins a set of points with 

dashed lines. The number and 
lengths of dashes and gaps are 
under program control, thus en- 
abling the user to draw any type 
of dashed lines. 

c. SPSYMB: See Section 3.4.4. 

d. PLTPLT: See Section 3.4.3. 



B. Geometric Figures : 

a. CIRCLE: Enables the user to plot a 

circle or a circular arc with a 
specified centre and radius. 

b. ELLIPS: Plots an ellipse or an 

elliptical arc with specified 
radii, centre and inclination. 

c. ARC: Draws a circular arc, or a 

circle, given the co-ordinates of 
the centre and the beginning and 
"end" of the arc. 

d. RECT: Plots a rectangle of speci- 

fied size with specified inclin- 
ation . 
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C. Grids and Axes : 

a. GRID: Plots a linear, semi-log or 

log-log grid of any size, anywhere 
on the ^plotting area with any re- 
quired density. 

b. AXES: Plots X-Y axes perpendicu- 

lar to each other anywhere on 
the plotting area, crossing at a 
specified point, with linear, 
semi-log or log-log tick marks 
which could be of any specified 
lengths . 

c. AXS: Plots a line, of specified 

length, at any angle with equally 
spaced tick marks of any speci- 
fied lengths and properly cen- 
tred annotation of any specified 
size on a specified side of the 
line. 

d. AXS 2 : Plots a line, between two 

specified points, with equally 
spaced tick marks of any speci- 
fied lengths and properly cen- 
tred annotation of any specified 
size on a specified side of the 
line. 



D. Results of Calculations : 

a. NMBR: Plots REAL, INTEGER, 

LOGICAL or COMPLEX number (s) of 
any length in a specified FORTRAN 
FORMAT at a specified location 
with any specified height and 
inclination. 



E. Miscellaneous : 

a. ROTPLT: Facilitates plotting with 

rotated X-Y axes. 

b. MIRPLT: Draws the mirror image of 

a complex plot by placing an im- 
aginary mirror at any desired 
angle. 



3.4.7 Basic Plot Routines : 

The X-Y co-ordinates of all the 
points given to the routines discussed so 
far are in the user's units defined by 
the parameters of the subroutines AREA 
and SETPLT. These parameters are re- 
ferenced between the subroutines through 
a labelled COMMON P611G. All these plot- 
ting routines directly or indirectly in- 
voke the plot routines which handle X-Y 
points in one thousandths of an inch. 
The formulae which define the relations 
between X,Y in user units and IX,IY in 
3ne thousandths of an inch, in terms of 
/arious parameters in the labelled 
COMMON P611G, are: 



IX=(X-XMIN)/XSCALE*XINCH 
IY=(Y-YMIN)/YSCALE*YINCH 

X=rX/XINCH*XSCALE+XMIN 
Y= I Y/YINCH+YSCALE+YMIN 

where , 

XMINl - Minimum X and Y values in 
YMINj user units. 

XSCALE*l- Lengths of the plotting area 
YSCALEJ in X and Y directions in the 
user units. 

XlNCHl - Physical dimensions, is one 
YINCHj thousandths of an inch, iri X 
and Y directions. 

The general user seldom has a need to 
call these basic routines which are: 

a. PLTRTN: Moves the plotting pencil in 

up or down position with X-Y co-or- 
dinates specified in one thousandths 
of an inch. 

b. LCTRTN: Returns the co-ordinates, in 

one thousandths of an inch, of the 
current position of the plotting 
pencil . 

c. SFTRTN: Enables the user to shift, by 

a specified amount in one thousandths 
of an inch, all future plotting with 
respect to plotting done prior to the 
call to SFTRTN. 



4.0 Backbone of the System : 

All the routines discussed so far are 
device independent in the sense that the 
user may call them while plotting on any of 
the plotting devices - pseudo or real. Now 
an attempt will be made to describe how 
this device independence is achieved. 

There are four plotting functions 
which depend upon the nature of the plot- 
ting device, namely, 

a. physical plotting, 

b. location of the current position of the 
plotting pencil , 

c. termination of a plot, and 

d. shifting plots. 

The above functions are achieved by the 
subroutines PLTRTN, LCTRTN, ENDPLT and 
SFTRTN respectively. (To be precise, 
LCTRTN, ENDPLT and SFTRTN are entry points 
of PLTRTN.) Though function b listed 
above is not strictly device dependent, it 
is convenient and straightforward to 
associate it with physical plotting. 
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All the plotting routines which pro- 
duce plots call the subroutine PLTRTN dir- 
ectly or indirectly. Thus, the device in- 
dependence of the entire plotting package, 
as far as a choice of the plotting device 
is concerned, is based on the availability 
of different versions of the subroutine 
PLTRTN (with entry points LCTRTN, ENDPLT 
and SFTRTN) . 



4 . 1 Selection of a Device t 

A typical FORTRAN job in IBM System 
370 goes through the following distinct 
processes (see Ref. 4 and 5): 

a. Compilation, 

b. Linkage editor (or loader) and 

c. Execution. 

During the first process, the FORTRAN 
program (including the supplied subpro- 
grams) is translated into machine language 
to produce what is called an "object 
module". The calls made to the non- 
supplied subprograms are generated as ex- 
ternal references and are marked as unre- 
solved. 

During the linkage editor (or loader) 
process, all the external references are 
resolved in the following manner: 

a. All the object modules given as input 
to the linkage editor (or loader) are 
merged with the object module produced 
as a result of the compilation. 

b. All the remaining external references 
are resolved from the supplied libra- 
ries (specified through the SYSLIB 
data definition card) . These libra- 
ries are searched in the given order. 
It is important to note that the 
search for a particular routine ter- 
minates as soon as it is found. 

The module produced as a result of the 
second process is referred to as a "load 
module". During the last process, the 
load module is executed. 

Several versions of the subroutine 
PLTRTN (with entry points LCTRTN, ENDPLT 
and SFTRTN) have been written and stored 
in various different libraries. The pro- 
per version of PLTRTN is obtained by pro- 
viding suitable concatenation of the data 
sets for library search through the SYSLIB 
DD card for the linkage editor or the 
loader. This process can be conveniently 
automated for a common user by creating 
appropriate JCL procedures and storing 
them in the system procedure library 
(SYS1.PROCLIB) . 

It follows from the above discussion 
that the addition of a new plotting device. 



whether it be a real or pseudo, is also a 
simple matter of writing the subroutine 
PLTRTN (with entry points LCTRTN, ENDPLT 
and SFTRTN) for that particular device. 
There need not be any change whatsoever 
either to the user's programs or to the 
rest of the plotting package. 

Testing of a new plotting device is a 
very straightforward process since the 
system programmer can supply the PLTRTN as 
a source program - either FORTRAN or 
assembler , as an object module or as the 
load module. This testing can proceed con- 
currently with the other user programs us- 
ing other plotting devices without any 
interaction. 



4 . 2 Various Plotting Devices t 

It would be beyond the length of this 
paper to describe all the plotting devices 
that have been and are connected to the 
system. Therefore, only a few major ones 
will be discussed in this section. 

In the off-line mode of plotting, the 
X-Y points are stored as ICNT, IX and IY 
where ICNT is the control byte and IX, IY 
are the coordinates of a point in one 
thousandths of ah inch. The meanings of 
the control byte are: 

ICNT ACTION 

0 move plotting pencil to (IX, IY) in 
up position 

1 lift plotting pencil up, then 
move it to (IX, IY) and put it down; 
i.e. plot the point (IX, IY) 

2 draw a straight line from previous 
plotting pencil position to (IX, IY) 
and leave the plotting pencil down 

3 do as ICNT=2 above and lift the 
plotting pencil ud 

4 call ENDPLT; ignore (IX, IY) 

5 call SFTRTN with (IX, IY) as the 
arguments. 

The plotting package arranges ICNT such 
that 0 s ICNT s 5 . 



4.2.1 Off-Line Plotting Through 
MASSPLOT : 

PLTRTN for this version is stored 
in the library UNB1- PLOT -MASSPLOT. This 
version stores the plots produced by the 
users in a predefined and formatted direct 
access data set. The user is never re- 
quired to create the data set. The plot- 
ting package stores several plot frames 
for several users in this data set which 
consists of 255 blocks each containing 
748 bytes. The first block is the direc- 
tory block and the remaining are plot-data 
blocks. 

An entry is created in the directory 
for each plot frame stored. The format of 
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this directory entry is given in Fig. 6. 
The first byte contains the plot frame 
number starting from 11 and the last byte 
contains the block number where the plot- 
data for this plot frame begins. Plot 
code consists of 11 packed decimal digits 
giving the time of the day the plot frame 
was created. The first byte of the plot 
code is made hexadecimal FF (=255iq ) when 
the plot is deleted. 



The control program gives the plotter oper- 
ator various options such as producing the 
plot frame on the specified device, skip- 
ping the plot frame, deleting the plot 
frame, forcing the plot frame on any of the 
available real devices, clearing the entire 
data set, etc. 



4.2.2 Off-Line Plotting Through 
USERPLOT: 



Byte number 

(* 1 2 3 4 5 6 


7 8 9 10 


1 1 1 1 1 1 






~- 1 


1 ^~ 

Plot 

frame 

number 


Plot Code 


Reserved | 

Block number 
where plot 
begins 



Fig. 6 Directory Entry for MASSPLOT 



The format of the data structure used 
for the plot-data blocks is given in Fig. 
7. 



First Plot-Data Block: 



Byte number 
12 3 4 5 
I I I I 1~~ 



c 



10 



I I I I 



Device 

1- 611 

2- 1627 

3- 5100 



XINCH 



YINCH 



N£80 



CByte number 
M+l M+2 M+3 M+4 M+5 MH> M+7 M+8 M+9 



ICNT 



IX 



IY 



If N*0, no user identification label 

if 0<N£80, bytes 11 to 10+N contain the user 

identification label 
M=10+9xf(N/9) where T gives the ceiling 
Remaining bytes contain 9 byte groups of 

ICNT, IX, IY 



Subsequent Plot-Data Blocks : 

C- Byte number 
12 3 4 5 6 7 



I 1 I I I I I 



8 

o 



10 



♦ * , 

Plot ICNT 

frame IX IY 

number 

Remaining bytes contain 9 byte groups of 
ICHT,IX,IY 



Fig. 7 Plot-Data Blocks for MASSPLOT 



PLTRTN for this version is stored 
in the data set UNBl • PLOT • USERPLOT . It 
stores plots for a particular user in a 
sequential data set created by him. Since 
the data set can contain the plot frames 
for only one user, no directory is required. 

Each record is 126 bytes long and has 
the format shown in Fig. 8. The only re- 
striction on the data set is that the re- 
cord length must be 126 bytes. It may be 
blocked and may be created on any type of 
available auxiliary storage device (e.g. a 
disk or a tape) . 



First Record of a Plot Frame: 



C- Byte number 
12 3 4 



126 



Fra m e I Plot No. 
Separator (01 to 99) 
(F0 )V EBCDIC 

At> Device N<80 

1- 611 

2- 1627 

3- 5100 

If N-0, remaining bytes are ignored 
if 0<Ns80, bytes 6 to 5+N contain the user 
identification label 

Subsequent Records : 

C- Byte number 
123456789 



till i 



ICNT 



IX 



IY 



Remaining bytes contain 9 byte groups of 
ICNT, IX, I Y 



Figure 8 Data Structure for USERPLOT 



The same data structure is used for 
HASPPLOT - in fact, USERPLOT has been 
created by a very simple and straightfor- 
ward simulation of HASP local modifications. 
That is precisely the reason why the record 
length is so small. 
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As in the MASSPLOT version, the con- 
trol program used to retrieve the plots 
stored in USERPLOT system gives the plot- 
ter operator various options such as pro- 
ducing the plot frame on the specified de- 
vice, skipping the plot frame, forcing the 
plot frame on any of the available real 
plotting devices, etc. 



4.2.3 On-Line Mode (REALPLOT) : 

PLTRTN for this version is stored 
in the library UNB1 • PLOT • REALPLOT . This is 
the version that actually produces plots. 

As created at the author's installa- 
tion, this version gives access to the 
following three devices: 

a. 611 storage oscilloscope, 

b. 1627 X-Y plotter and 

c. 5100 electrostatic plotter 

The choice is governed through a full word 
integer in the labelled COMMON P611G. 

Basically, the IX, IY values passed to 
PLTRTN or SFTRTN are scaled to the proper 
units for the selected device and these 
are then passed on to the driver routines 
(see Section 5.0} which do the plotting. 
ENDPLT simply invokes the appropriate 
driver routine. Since the information 
required for LCTRTN is already in PLTRTN 
itself, no corresponding driver routines 
are required. 

The control programs for MASSPLOT 
and USERPLOT invoke this version of PLTRTN 
to produce the actual plots for the user. 

It is perhaps quite obvious by now 
that one could create versions of PLTRTN 
which plot only on the 611 or the 162 7 or 
the 5100 and store them in three different 
libraries. Then, the device selection 
can be achieved even in the on-line mode 
of plotting by a simple JCL change. 



4.2.4 Printer Plots : 

PLTRTN for this version is 
stored in the library UNBl -FORTLIB. Since 
this is the default library searched for 
resolving external reference for the link- 
age editor or loader process in all the 
FORTRAN procedures at the author's instal- 
lation, the printer plots can be produced 
without any special JCL cards. 

The version currently used in produc- 
tion builds the plot in a two dimensioned 
array which is printed when ENDPLT is 
called. This forces the arbitrary limits 
on the size of the plot that can be pro- 
duced on the printer. 



Several other versions have been writ- 
ten and successfully used which remove the 
arbitrary restriction on size and which 
use various different algorithms. 



5.0 Device Dependent Routines : 

Though the device dependence may be 
handled in PLTRTN itself, it is often con- 
venient, and indeed sometimes desirable in 
order to access several real plotting de- 
vices during the execution of a program, to 
localize the device dependency in the 
driver routines. These driver routines may 
then be called by PLTRTN. 



5.1 Drivers : 

The driver routines are: 

Pxxxx - movement of plotting pencil 

Sxxxx - shifting of plots 

Exxxx - termination of plots 

where, xxxx is the device. PLTRTN calls 
Pxxxx, SFTRTN calls Sxxxx and ENDPLT. calls 
Exxxx. 

Thus, the routines for the 611 storage 
oscilloscope aw: 

CALL P61K IC, IX,1Y) 
CALL S611UX,IY) 
CALL E611 

where IX and IY are in 611 absolute units, 
i.e. between 0 and 32742. 

The routines that drive 1627 and 5100 

are: 

CALL P1627I IC f IX » I Y ) 
CALL S1627UX.IY) 
CALL E1627 

CALL P5100UC,IX,IY) 
CALL SSIOOUX.IY) 
CALL E5100 

where IX and IY are in one hundredths of 
an inch. 



5.1.1 Tektronix 611 Storage 
Oscilloscope : 

This is a storage cathode ray 
tube with a physical picture area of 6.37 5 
inches in the X direction and B.125 inches 
in the Y direction. The full scale deflec- 
tion is produced by applying one volt to 
the X and/or Y direction. Since the 611 is 
controlled by the IBM 370/158 computer 
through the 1827 data control unit, the 
number 32742 generates one volt and creates 
the full scale deflection. 
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The algorithm used to generate 
straight lines is an extension of symmetri- 
cal digital differential analyzer (ESDDA) . 

Assume that the beam can store its 
path if up to MAXINC increment is given 
in the X and /or Y direction. MAXINC may 
be typically 64 or even higher. Then the 
beam can be moved in MAXINC different ways 
in any octant. See Fig. 9 which depicts 
various possible movements for different 
values of MAXINC. By choosing the move- 
ment that is closest to the desired 
straight line, one can obtain a fairly 
good approximation. 




<c) MAXINC-4 

Fig. 9: Possible Movements for Various 
Values of MAXINC 



The flowchart of the algorithm, 
based on this idea to draw a straight 
line from (xo,YO) to (x,Y) is given in 
Fig. 10. The following notation is 
used in Fig. 10: 

| d | ^absolute value of d 

REM (d , x) ^Remainder of d*x 

SIGN (design of d 

The reader might recognize it as the gen- 
eralized symmetrical digital differential 
analyzer, the later being a special case 
when MAXINC is one. 

5.1.2 Calcomp 1627 X-Y Plotter : 

As is perhaps well known, the 
Calcomp 162 7 plotter works in steps of 



one-hundredths of an inch and as such the 
drivers work in that unit. 

The algorithm used is the one given in 
Ref. 6. Concise mathematical definition of 
this algorithm appears in Appendix C. 



5.1.3 Gould 5100 Electrostatic 
Plotter : 

The algorithm developed for this 
device is the subject matter of another 
technical report (see Ref. 7). 



5.2 Odds and Ends : 

There are certain other types of de- 
vice dependencies introduced by the pecu- 
liar nature of the plotting device. For 
example, subroutines PRNTCH, CHLNIN and 
CHRPRT are required for line printers. It 
would be a great inconvenience, and indeed 
a great nuisance, if the user were required 
to insert and remove certain statements 
from his program depending upon the plot- 
ting device. Appropriate dummy subroutines, 
requiring sometimes as few as two bytes, 
are written and placed in proper libraries 
to take care of such situations. 



6.0 Plot Accounting : 

The method used for charging the plots 
produced is given in this section. This 
method, which was designed for the 611 
storage scope and the 1627 X-Y plotter, may 
undergo some modifications for the 5100 
electrostatic plotter. 

The effective plot time (EPT) is cal- 
culated according to the following formula: 

EPT=A*PCLS+B*PINCH milliseconds 

where , 

PCLS - Number of calls to PLTRTN 
PINCH - Plot inches 

A - milliseconds required to lift 

the plotting pencil up or put 
it down 

B - milliseconds required to move 
the plotting pencil one inch 

PCLS is a good approximation of up and 
down movements of the plotting pencil. 

PINCH is calculated according to the 
following formula: 

PINCH = I ( |X-XO[ + |Y-Y0|) 

where (XO,YO) and (X,Y) are the starting 
and end points of each plotting pencil 
movement and the summation is over all 
motions used to produce a plot frame . 
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SWITCH*-* TRUE • 
INCRY*ttAXINC*SIGN (DY) 

length*- 1 dy | 
incrx+maxinc*dx/ 1 dy , 
remndrh rem (maxinc*dx , i dy i ) 
valcor*-|dy| 



INTEGER X,Y,XO,YO,DX,DY 

INCRX , INCRY , LENGTH , 
VALCOR, COUNT , REKNDR , 
MAXINC/64/ 

LOGICAL*! SWITCH 



No 



1 



SWITCH*- -FALSE • 
INCRX*-MAXLNC*SIGN (DX) 
LENGTH*- 1 DX | 
INCRY*^MAXINC*DY/ | DX | 
REMNDR* I REM (MAXINC*DY , | DX | ) | 
VALCOR*- 1 DX | 



COUNT *- VALCOR/ 2 



< 



Hove to (XO.YO) 



> 



XO +■ XO + INCRX 
YO <- YO + INC BY 
COUNT *- COUNT + REMNDR 



COUNT *- COUNT - VALCOR 




XO*-X0-l 



XO*-XO+l 



YO*-YO-l 



YO-YC+1 



— 



LENGTH *- LENGTH - MAX INC 



Yes 



0 — < 




The portion in this dotted 
box can be programmed with 
one instruction in 360 
assembler* 



> 



Fig. 10: Algorithm for generating a straight line between (XO.YO) and (X,Y) 
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The values of constants A and B depend 
upon the plotting device. The current 
values are: 



Device 


A 


B 


611 
1627 


.1 

140 


17 
500 



The plotting charge (PC) is given by 
the following formula: 

PC «* C*EPT4-D units 

where C and D are constants which are cur- 
rently set at 0.00025 and 40 respectively. 



7.0 Conclusions : 

A powerful computer plotting system 
has been presented. The system has been 
in use for a number of years and has 
proved to be satisfactory to several users. 

Several new plotting devices have 
been added since the inception of the sys- 
tem. These additions have been completely 
transparent to the users. 

The aspects that are dependent on the 
plotting device and those, that are indepen- 
dent of the plotting device have been 
separated and identified. Based on this, 
the interface for adding a new plotting 
device has been formalized and defined. 
The implementation of this interface has 
been straightforward and has worked very 
satisfactorily from the system program- 
mer's point of view. 

The package is open ended in the 
sense that it can be expanded both in the 
devide independent section as well as de- 
vice dependent section. For example, 
addition of a contouring package or three 
dimensional plotting once implemented will 
work for all the plotting devices. 

As an example of a new pseudo plot- 
ting device, one may consider the imple- 
mentation of two dimensional transforms 
(see Ref, 8) to obtain the effects of 
translation, rotation, scaling, etc. 
Further, one could incorporate the 
windowing, clipping, etc. 
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Appendix A : Summary of Calling Sequences 
of the Routines 



In the following list, the optional argu- 
ments are. enclosed in square brackets. 
The routines marked with an- asterisk (*) 
are the drivers and as such are device 
dependent; all other routines are device 
independent. 
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CALL ARC<XC,YC,XS,YS[,XL,YL] i 
CALL AREA(XLNGTH.YLNGTH) 

CALL AXES(XMN,YMN,XMX,YMX,XCROSS,YCROSS, 

XDt YD( tTICKLN] > 
CALL AXSU.Y, I ARRAY .SIZE, THETA, VALUE, 

VALINC ,NDD [ i HIGHTL»HIGHTNf TICKLN] ) 
CALL AXS2(X1,Y1,X2,Y2, 1 ARRAY, VALUE, 

VALINC, NDD( ,HIGHTL,HIGHTN , TICKLN] ) 
CALL CHARS (X,Y, I ARRAY [, THE T A , HEIGHT t 

XRET , YRET J ) 
CALL CHLNIN ( XCHARS , YLNS ) 
CALL CHRPRT(X i Y, I ARRAY t , THETA , HEIGHT, 

XRET, YRET] ) 
CALL CIRCLE(XC,YC,RADIUS [ , FROMTH , TOTH] J 
CALL DEVICE(ID) 

CALL DSHLNS(XX,YY[,DSHGPL,IDG, ILIM, 
ISTART,IINCR]) 

CALL ELLI PS ( XC » YC ,RMAJOR [ , R MI NOR , THETA i 
FROMTH, TOTH] ) 

CALL ENLPLT 
♦CALL E1627 
*CALL E5100 
♦ CALL E611 

CALL GR I D ( XMN , YMN , XMX , YMX , XD , YD ) 

CALL LCTRTN(IX,IY) 

CALL LNSPLT(XX,YY[,ILIM,ISTART,IINCR] ) 
CALL LOCATE(X,Y) 

CALL HIRPLTUC,X,Y,THETA[ ,XHIR,YMIR] ) 
CALL NMBR (X,Y, NUMBER, FORMAT [, THETA, HEIGHT, 

ILIM, I START, II NCR, XRET, YRET] ) 
CALL NOWPLT( IC,X,YJ 
CALL PLOTID( I ARRAY ] 

CALL PLTPLT(ICC,XX,YY[,ILIM,ISTART,IINCR] ) 

CALL PLTRTN(IC,IX,IY) 

CALL PRNTCH(ICH) 
♦CALL P1627( ICIXtlY) 
♦CALL P510O( IC,IX,IY) 
♦CALL P611UC, IX, IY} 

CALL RECT 1 X»Y t XL [ ,YL, THETA] ) 

CALL ROTPLTt IC » XROT i YROT ,THETA[,X,Y] ) 

CALL SETPLT(XMN,YMN,XMX,YMX) 

CALL SFTRTNt IX, IY) 

CALL SHIFT(X,Y) 

CALL SPSYM3 tX,Y, I [, THETA, HEIGHT, XRET, YRET] ) 
♦CALL S1627( IX, IY) 
♦CALL S5100(IX,IY) 
♦CALL S611(IX,IY) 

CALL THlCMXrY, I ARRAY, THETA, HEIGHT, IDX,IDY> 



Appendix B : Libraries and Plotting Routines 

Notes: 1. P611G referred to in the following 
tables is a labelled COMMON block.- 
2. F and A in .the column "written in B 
of the following tables stand for 
FORTRAN IV and 370 Assembler 
respectively. 



A. Library: UNB1 -FORTLIB 



Routines 


Entry 
Points 


Memory 
(bytes) 


Other 
Routines 
Called 


Written 
In 


ARC 




1058 


VARARG 

PLTRTN 

ATAN2 

SQRT 

COS 

sra 

P611G 


F 


AREA 




334 


P611G 


F 


AXS 


AXS2 


3892 . 


VARARG 

CORE 

IBCOH0 

CHRPRT 

PLTRTN 

AT AN 2 

SQRT 

COS 

SIN 

P611G 


F 


CHARS 


SPSYMB 


2864 


PLTRTN 
COS 
SIN 
P611G 


A 


CIRCLE 




946 


VARARG 

PLTRTN 

COS 

SIN 

P611G 


F 


DEVICE 




2 




A 


DSHXNS 


- 


1698 


VARARG 

LOCATE 

PLTRTN 

SQRT 

P611C 


F 


ELLIPS 




1188 


VARARG 

PLTRTN 

COS 

SIN 

P611G 


F 


GRID 


AXES 


2720 


FIXPRtf 
FlXPItf 
NOWPLT 
LOCATE 
VARARG 
ALOGIO 
P611C 


F 


LOCATE 




424 


LCTRTN 
P611G 


F 


MIRPLT 




732 


VARARG 

PLTRTN 

COS 

SIN 

P611G 


F 
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Library: UNB1 -TORTLIB (continued) 



C. Library: UKB1 ■ PLOT-USERPLOT 



Routines 


Entry 
PointB 


Memory 
(bytes) 


Other 
Routines 
Called 


Written 
In 


NMBR 




1188 


tff A ft a n *^ 

VARARG 

CORE 

IBCOHff 

CHRPRT 

P611G 


F 


NOWPLT 




458 


PLTRTN 
P611G 


F 


PLOTID 




2 




A 


PLTPLT 


LNSPLT 


444 


PLTRTN 
P611G 


A 


PLTRTN 


CHRPRT 
CHLNIN 
ENDPLT 
LCTRTN 
PRNTCH 
SFTRTN 


4082 


IBCOM0 
VARARG 
COS 
SIN 

P1403G 
P611G 


F 


RECT 




1266 


VARARG 

PLTRTN 

COS 

SIN 

P611G 


F 


ROTPLT 


- 


708 


VARARG 

PLTRTN 

COS 

SIN 

P611G 


F 


SETPLT 




394 


P611G 


F 


SHIFT 




428 


SFTRTN 
P611G 


F 


THICK 




1898 


CHARS 
COS 
SIN 
P611G 


F 



Library: UNBl-PLOT'MASSPLOT 



Routine 


Entry 
Points 


Memory 
(bytes) 


Other 
Routines 
Called 


Written 
In 


CHRPRT 




12 


CHARS 


A 


DEVICE 




124 


P611G 


A 


PLTRTN 


ENDPLT 
LCTRTN 
PLOTID 
SFTRTN 


2840 


P611G 

Acctng. 

modules 


A 


PRNTCH 


CHLNIN 


2 




A 



Routine 


Entry 
Points 


Memory 
(bytes) 


Other 
Routines 
Called 


Written 
In 


CHRPRT 




12 


CHARS 


A 


DEVICE 


- 


124 


P611G 


A 


PLTRTN 


ENDPLT 
LCTRTN 
PLOTID 
SFTRTN 


2532 


WTL 

IBCOM# 

P611G 

Acctng. 

modules 


F 


PRNTCH 


CHLNIN 


2 




A 



D. Library: UNB 1 * PLOT • REALPLOT 



Routine 


Entry 
Points 


Memory 
(bytes) 


Other 
Routines 
Called 


Written 
In 


CHRPRT 




12 


CHARS 


A 


DEVICE 




124 


P611G 


A 


PLOTID 




2 




A 


PLTRTN 


ENDPLT 
LCTRTN 
SFTRTN 


1518 


P611 
P1627 
P5100 
E611 

E1627 

E5100 

S611 

S1627 

S5100 

Acctng. 

modules 


F 


PRNTCH 


CHLNIN 


2 




A 


P1627 


E1627 
S1627 


1052 




A 


P5100 


E5100 
S5100 


7728 


SORT 


A 


P611 


E611 
S611 


1852 




A 



Api>endix C : Algorithm Used for the 1627 

The following is a concise mathematical definition 
of the algorithm by Bresenham (Ref . 6) to draw a 
straight line from (XO,YO) to (X,Y) on the 1627 
X-Y plotter. 



If the possible plotter movements are: 



j + ve Y 




+ve X 



then the algorithm is: 

1. Compute DX and DY as 

DX = X-XO and DY ~ Y-YO 
Determine DA and DB as 

DA = Maximum of |DX| and |dy| 
DB = Minimum of | DX | and | DY | 
Calculate DEL as 2*DB-DA 
Compute Ml and M2 as 



2. 



3. 
4. 



5. 



6. 



Ml = 7 


if 


DX 




DY 


and 


DX < 0 


= 3 


if 


DX 




DY 


and 


DX a o 


= 1 


if 


DX 


< 


DY 


and 


DX £ 0 


= 5 


if 


DX 


< 


DY 


and 


DX < 0 



M2 = 



2 if 
4 if 
6 if 
8 if 



DX 
DX 
DX 
DX 



£ 0 and DY £ 0 

£ 0 and DY <:0 

< 0 and DY < 0 

< 0 and DY £ 0 



If DEL2O , execute the plotter move- 
ment M2 and modify DEL as: 

DEL-DEL+2*DB-2*DA 

If DEL<0 , execute the plotter move- 
ment Ml and modify DEL as: 

DEL=DEL+2*DB 

Repeat Step 5 until you reach 
<X,Y). 
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